Absence of the discontinuous transition in the one-dimensional triplet creation model 
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Although Hinrichsen in his unpublished work theoretically rebutted the possibility of the discon- 
tinuous transition in one-dimensional nonequilibrium systems unless there are additional conserva- 
tion laws, long-range interactions, macroscopic currents, or special boundary conditions, we have 
recently observed the resurrection of the claim that the triplet creation model (TC) introduced by 
Dickman and Tome [Phys. Rev. E 44, 4833 (1991)] would show the discontinuous transition. By 
extensive simulations, however, we find that the one-dimensional TC does belong to the directed 
percolation universality class even for larger diffusion constant than the suggested tricritical point 
in the literature. Furthermore, we find that the phase boundary is well described by the crossover 
from the mean field to the directed percolation, which supports the claim that the one-dimensional 
TC does not exhibit a discontinuous transition. 

PACS numbers: 05.70.Ln,64.60.De,05.70.Fh 
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I. INTRODUCTION 

Although the field theory for the tricritical phenomena 
in reaction-diffusion systems was developed more than 
two decades ago [11, 0] (see also Ref . ) , not many nu- 
merical studies have followed [1, Q . One apparent reason 
is the numerical difficulty, but it can be soon overcome by 
the increasing computing power. More seriously, it was 
strongly argued that no discontinuous transition is pos- 
sible in one dimension once there are no additional con- 
servation laws, long-range interactions, macroscopic cur- 
rents, or special boundary conditions [6|, which rebutted 
the observed discontinuous transition in one-dimensional 
triplet creation model (TC) by Dickman and Tome Q. 
In view of the fact that large portion of the studies on 
the absorbing phase transitions (for a review, see, e.g. 
Refs. is focused on the systems in one dimen- 

sion, this theory presumably has kept researchers from 
being into the tricritical phenomena. 

Recently, however, numerical studies in favor of the 
original claim by Dickman and Tome have been re- 
ported [n], [12, [iJl . If this claim rather than the theory 
in Ref. '6!| turns out to be right, we would observe an 
avalanche of studies on the tricritical phenomena. Un- 
fortunately, however, no theoretical argument regarding 
the mechanism to stabilize a domain in one dimension has 
been suggested as yet. Moreover, the tricritical point of 
the diffusion rate reported in Ref. [ij, [13] is too large 
to reject the opinion that the system will eventually 
crossover to the directed percolation (DP) universality 
class after a long transient time. Recent study by Ferreira 
and Fontanari |I4|] using n-site approximation alluded to 
the crossover rather than the tricritical behavior, though 
they did not strongly put forward such a scenario because 
of the computational limitation of their method. Actu- 
ally, Hinrichsen [Ql numerically showed that the simula- 
tion time in Ref. [7| was too short to sec the correct scal- 
ing behavior. Interestingly, Cardozo and Fontanari (l^ 
also refuted the value of the tricritical point originally 
suggested. Hence, if the argument by Hinrichsen Q is 



right, it is very probable that more extensive simulations 
than those in Ref. [l3| would revive the history. 

Indeed, we found the DP scaling over the parameter 
range where the tricritical point was located by Cardozo 
and Fontanari p^ . This paper is for providing numer- 
ical evidences to support the theory suggested by Hin- 
richsen 0. 

The rest of this paper is organized as follows: Sec- 
tion |TT] introduces the d-dimensional TC and explains the 
algorithm implemented for numerical simulations. The 
numerical results showing the DP scaling behavior for 
larger diffusion rate than the previously reported tricriti- 
cal point will be presented in Sec. lIIII In Sec.|TVl we will 
argue that there is a crossover rather than a tricriticality 
in one dimension by studying the behavior of the phase 
boundary. Section IVl summarizes the work. 



II. TRIPLET CREATION MODEL 

The TC is an interacting hard core particles system on 
a d-dimensional hypercubic lattice with three processes, 
hopping (with rate D), spontaneous annihilation (with 
rate 7), and creation by a triplet (with rates s) [7]. By 
hard core is meant that no two particles can occupy the 
same site at the same time. By a suitable time rescaling, 
we can set D + 'y + s — I without loss of generality. It 
is also convenient to introduce the annihilation probabil- 
ity p such that 7 = (1 — D)p and s = (1 — D){1 — p). 
The detailed dynamics is to be explained in terms of the 
algorithm used for simulations. 

At time t, Nt particles are distributed on a d- 
dimensional hypercube of size L"* [Nt is a random vari- 
able). Each site is represented by a lattice vector m ~ 
(toi, . . . , md) {0 < rrii < L — 1). The unit vector along 
direction i is denoted by (i = 1 , . . . , d) . In all simula- 
tions in this paper, a fully occupied initial condition and 
periodic boundary conditions are assumed. 

The algorithm begins with a random selection of a par- 
ticle in the system. For convenience, let us refer to the 
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lattice vector of the site the selected particle resides as 
m. With probability D, hopping is attempted to a target 
site which is chosen randomly among 2d nearest neigh- 
bors of the site m. This hopping is successful only if 
the target site is empty (hard core exclusion), otherwise 
there is no configuration change. With probability l — D, 
either annihilation (with probability p) or creation (with 
probability 1 — p) will be attempted. When annihilation 
is decided, the selected particle will be irreversibly re- 
moved from the system. If the creation is to occur, one 
of the directions i {i — l,...,c?) is selected randomly. 
Two sites m-\-ei and m + 2ei are checked whether both 
sites are occupied or not. If both sites are also occupied, 
one of the two sites, m + Se; or m — e;, is chosen at 
random as a target site and a new particle is created at 
the target site provided it is empty. If any of the condi- 
tions for creation is not satisfied, nothing happens. After 
an attempt to change a configuration, time increases by 
l/Nt regardless of its success. The above procedure iter- 
ates until either the system reaches the absorbing state 
where no particle remains or time gets larger than the 
preassigned observation time. 

When d = 1, the above dynamics is exactly the same 
as that in Ref. [ll] with p = 1/(1-1- A). This is because 
the probability to find three occupied sites in a row does 
not depend on whether two sites m + Si and m + 2ei 
(in this work) or m + Bi and m — Bi (in Ref. [l2|) are 
examined. 

Although we introduced d-dimensional TC, discus- 
sions from Sec. IIIII on will be focused only on the one- 
dimensional model. Thus, the dimensionality of the 
model will not be mentioned explicitly in the following. 

In the simulation, we measure the particle density 
p{t) = {Nt)/L, where (. . .) means the average over inde- 
pendent realizations, and the (survival) probability that 
there is a particle in the system at time t. The measure- 
ment of the survival probability is mainly for making sure 
that the system is large enough not to be affected by the 
finite size effect up to the observation time. 



III. CRITICAL DENSITY DECAY 

This section studies the critical behavior of the TC 
with D = 0.95 and D = 0.98. Since the question in 
this section is whether the TC shows the discontinuous 
transition or the continuous transition governed by the 
DP fixed point, observing the DP scaling behavior of a 
single quantity is enough for our purpose. Anticipating 
the conclusion, we only study how the density decays 
near criticality. 

Figure [1] depicts p{t)t^ as a function of t in semi- 
logarithmic scales, where 5 = 0.1595 is the critical den- 
sity decay exponent of the DP class taken from Ref. [l^ . 
The system size in the simulations for D — 0.95 (0.98) is 
L — 2^^ (2^^). The number of independent runs for each 
data set varies from 16 {D = 0.95 and p = 0.089 89) to 
100 {D = 0.98 and p = 0.094 226). The system evolves 
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FIG. 1: (color online) Plots of p{t)t" vs t with 5 = 0.1595 
(critical exponent of the DP class) for D = 0.95 (lower three 
curves) and D — 0.98 (upper three curves) in semi-logarithmic 
scales. The values of p for D = 0.95 are 0.089 89, 0.089 895, 
and 0.0899 from top to bottom. The values of p for D = 0.98 
are 0.094 224, 0.094 226, and 0.094 228 from top to bottom. 



up to i = 10^ at the longest and no sample falls into 
the absorbing state during simulations. For D = 0.95, 
the curve corresponding to p = 0.089 89 (0.0899) veers 
up (down), which indicates that the system is in the ac- 
tive (absorbing) phase. At p = 0.089 895, the curve is 
flat for more than two log-decades. Hence we conclude 
that the TC with D = 0.95 belongs to the DP class 
with critical point pc = 0.089 895(5), where the num- 
ber in parentheses indicates the error of the last digit. 
If we write the critical point using 7 = (1 — D)p, we 
get 7c = 0.004 4948(3) which should be compared with 
0.004 50(1) reported in Ref. [12]. One should note that 
the DP scaling is observable from t = 10^ which is the 
end point of the simulation for D = 0.95 in Ref. [l3| 
(see Fig. 3(b) of Ref. [l^). Likewise, the simulation re- 
sults for D = 0.98 show the similar behavior as those for 
D = 0.95 (see upper three curves in Fig.[T]). The critical 
point for D = 0.98 is found to be pc = 0.094 226(2) or 
7c = 0.001 888 452(4). In Ref. the critical value 7c 
for D = 0.98 was reported as 0.001 886(2) and the simu- 
lation was terminated around i = 10^ from when the DP 
scaling is observable. 

To conclude this section, the TC \yp to D — 0.98 be- 
longs to the DP class and previous claim of the existence 
of the tricritical point below D = 0.98 is refuted. Our 
results also explain why Cardozo and Fontanari [l3| ob- 
served continuously varying exponents as well as the com- 
pact growth; the system was analyzed before the correct 
scaling behavior was detected. 



IV. CROSSOVER FROM THE MEAN FIELD 
TO THE DIRECTED PERCOLATION 

In Sec. mil we numerically confirmed that up to D = 
0.98 the TC does show continuous transition governed 
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by the DP fixed point. Although this refuted the pre- 
vious claim [ill . [l2i | that the transition nature changes 
at a certain D smaller than 0.98, the possibility of the 
discontinuous transition at 0.98 < I? < 1 is still open. 
To provide an evidence that the DP fixed point governs 
the critical behavior beyond D = 0.98, this section stud- 
ies the phase boundary near D = I with the focus on a 
possible crossover. 

The investigation of the phase boundary to settle a 
controversy is not without precedent. The present au- 
thor and his collaborator tried to resolve the controversy 
around the pair contact process with diffusion (for a re- 
view, see, e.g., Refs. [3, [13) by studying the phase 
boundary of the crossover models [ll|. Although this 
study could not elicit a full consensus, it certainly gives 
a hint about the system. So we think it is worth while 
to investigate the phase boundary of the TC. But, in 
Sec, mil we have presented the simulation results only for 
D = 0.95 and D = 0.98, which is certainly not enough 
to see the structure of the phase boundary. For a better 
bird's-eye view, we will include some other critical points 
for D < 0.95 even though we will not be presenting de- 
tails. 

Before directly analyzing the phase boundary, we will 
discuss some features of it which can be inferred without 
resorting to the nature of the transition. First, the tran- 
sition point Pc{D) is argued to be an increasing function 
of D. Second, Pc{D) is expected to approach to the tran- 
sition point of the mean field (MF) theory as D — > 1. By 
the MF theory in this paper is exclusively meant the one- 
site approximation (for a detail, one may consult Sec. 3 
of Ref. dl) 

^ = -pp+{l-p)p^l-p), (1) 

which exhibits a discontinuous transition with the tran- 
sition point po = 

To argue that Pc{D) increases with D, it is convenient 
to modify the TC by rescaling time r = (1 — D)t. Note 
that for D < 1, however close 13 is to 1, the steady state 
property of the modified model is identical to the origi- 
nal one. In this modified model, a single particle dies out 
with rate p and a triplet attempts to branch an offspring 
with rate 1—p. Obviously, these rates are not dependent 
on D. The diffusion rate is now D = D/(l — D) which 
is an increasing function of D. Since the annihilation oc- 
curs regardless of the environment of a to-be-annihilated 
particle, the diffusion cannot directly affect the particle 
number fluctuation due to the annihilation. On the other 
hand, the branching is influenced by the diffusion, since 
a triplet can be either broken or newly formed by the 
movement of particles. Hence, whether the diffusion en- 
hances or reduces the activity of the system for given p 
can be answered by understanding if the diffusion will in- 
crease or decrease the number of active triplets, namely, 
triplets with a vacant neighbor. 

Due to the hard core exclusion, a cluster of particles 
will lose one by the diffusion only at boundaries. How- 



ever, even if a cluster loses a single particle at a boundary, 
the number of active triplets does not decrease if the size 
of the cluster is larger than 3 (it can even increase the 
number of active triplets if the configuration change due 
to the diffusion is like 01101111 01110111, where 1 
stands for a particle and for a vacancy). Thus, the 
reduction of the number of active triplets by the diffu- 
sion occurs in very restricted situations. On the other 
hand, the diffusion can mediate the formation of (active) 
triplets in the bulk of region where none exists. To sum 
up, the diffusion tends to enhance the activity, which en- 
tails that the transition point should increase with D. 
Next, we will discuss the limiting value of Pc{D) as 
^ 1, which will be denoted by p. To find p, we 
start from arguing that the TC under the D ^ 1 limit is 
deeply related to the MF theory. This connection makes 
sense only when time is suitably rescaled as before. With 
this time scale, the diffusion rate D grows indefinitely as 
£) — > 1, but the annihilation and creation rates remain 
finite. From now on, we will call the limit D — > 1 the fast 
diffusion limit and the modified TC with rescaled time is 
always assumed when we are discussing the fast diffusion 
limit. 

A finite system under the fast diffusion limit can 
be interpreted as follows: Right after any reaction (ei- 
ther annihilation or creation), the system arrives at the 
steady state of the diffusion process in no time and re- 
mains there until another reaction occurs. In the steady 
state of the diffusion-only system, all possible configu- 
rations for given number of particles have equal proba- 
bility. Hence, the probability that a site is vacant and 
its three consecutive neighbors are occupied at time r 
is [L — n)(n)3/(L)4, where L is the system size and n 
is the number of particles in the system at time r, and 
(m)fc = m{m— 1) . . . (m — k + l). Of course, the probabil- 
ity that a site is occupied at time r is n/L. Since these 
probabilities do not depend on where the site is located, 
the probability distribution is fully specified by the num- 
ber of particles and the master equation of the TC under 
the fast diffusion limit is reduced to 

+ Cn-lPn-1 - {an + C„)P„, (2) 

where a„ = pn, c„ = (1 — p){L — n){n)3/{L — 1)3, and 
Pn , though the argument of it is not written explicitly, is 
the probability that there are n particles at time r. For 
convenience, we set Pl+i = P-i = 0. 

The meaning of a„ and c„ in Eq. ([2]) can be interpreted 
as follows: a„ means each particle (there are n particles) 
dies with rate p and c„ means that each vacant site (there 
are L — n vacant sites) becomes occupied with rate (1 — 
p)(n)y,/ [L — 1)3, which is the dynamic rule of the TC on 
a fully connected graph. Hence, the master equation of 
the TC on a fully connected graph is exactly Eq. If 
we take the thermodynamic limit {L — )■ c») to Eq. ([2]), 
the law of large numbers has the density follow Eq. ([T]), 
that is, the MF theory. Now it is clear why we used 
T in Eq. ([T]) as a time parameter. For convenience and 
because of an obvious reason, we will refer to the TC on 
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FIG. 2: (color online) Comparison of the numerical solutions 
of Eq. ^ with the simulation results of the TC with 1 — 
D = 10~*. As explained in the text, the time is rescaled as 
r = {l~D)t. The systems sizes are L — 2^ (two below curves, 
though indiscernible) and 2* (two above curves), p is fixed at 
i. For L = 2®, no difference is detectable between simulation 
and the MF solution. On the other hand, the system with 
i = 2* is distinct from the MF solution after r = 50 (MF 
solution is slightly above the simulation results). 



a fully connected graph with finite size as the MF model. 

Although the fast diffusion limit was taken to arrive at 
the MF model in the above discussion, one can observe 
the behavior of the MF model even for nonzero 1 — D 
once (1 - D)L'^ < 1 0. In Fig. ^ we compare the 
simulations of the TC for 1 — D = lO^"* to the (numerical) 
solutions of the MF model a,t p = |. The system sizes 
are L = 2^ and L = 2^. For the solutions of the MF 
model, we numerically integrate Eq. ([2]) with the initial 
condition P„ (r = 0) = 6nL here is the Kroncckcr delta 
symbol). The simulation results are not discernible from 
the behavior of the MF model for L = 2^ [(1 - D)L'^ w 
0.4], but clear deviation is observed when L = 2^. 

In the above discussion, the thermodynamic limit, 
when necessary, is preceded by the fast diffusion limit. 
However, what we are interested in is the behavior of the 
TC with infinite size under the fast diffusion limit. That 
is, the thermodynamic limit should be taken before the 
fast diffusion limit. When two limits are involved in the 
calculation, one should be careful about which limit is 
taken first. Fortunately, the order of two limits are irrel- 
evant in most cases. What is meant by most cases will 
become clear in due course. 

Our discussion commences with the limiting behavior 
of the TC in the active phase {p < p). Let the critical 
diffusion probability for given p be denoted by Dc{p). 
Since Pc{D) is a monotonous function, Dc{p), if exists, 
is uniquely determined. If 1 > Z) > Dq, where Dq is an 
arbitrary number strictly larger than Dc(p), the system is 
in the active phase and the correlation length is bounded 
for all D in this regime. If the system size is much larger 
than the bound of the correlation lengths in this regime, 
the behavior in the thermodynamic limit is observable for 
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FIG. 3: (color online) Comparison of the numerical solutions 
of Eq. IT} with the simulation results of the TC with 1 — D — 
10"^ and 1-D = 10"* at p = 0.094 226. As in Fig.0 r is the 
rescaled time. The systems size for simulations is L = 2^^. As 
1 — D decreases, the density approaches to the MF solution. 



any Dq < D < 1 even though the system is finite. Since 
we can always choose D such that (1 — D)L'^ <C 1 for 
given L, the TC under the fast diffusion limit preceded 
by the thermodynamic limit is identical to the MF theory 
at the same value of p. 

To support this argument. Fig. [3] compares the MF 
theory with simulations of the TC for D — 0.999 and 
D = 0.9999 at p = 0.094 226 which is the critical point 
for D — 0.98. Obviously, the density of the TC obtained 
from simulations approaches to the MF theory as 13 — * 1 . 

The above consideration reveals that p should not be 
larger than po, otherwise the steady state density for pq < 
p < p would decrease with D, which is contradictory to 
the role of the diffusion as a enhancer of the branching. 

li p > Po where the TC as well as the MF model is in 
the absorbing phase, we can arrive at the same conclusion 
as above, that is, the behavior of the TC under the fast 
diffusion limit is not affected by whether before or after 
the thermodynamic limit is taken. 

To complete our discussion, we first have to figure out 
if p can be strictly smaller than pq. Actually, this possi- 
bility was suggested by Fiore and de Oliveira [ll| by ex- 
trapolating the phase boundary which was obtained nu- 
merically. Note that the discontinuity of a phase bound- 
ary per se is not an unrealistic conclusion. One can even 
find another report which, though in a different context, 
shows the discontinuity of the phase boundary [l^ . 

However, this scenario does not seem plausible. To 
demonstrate why it is not likely, let us think about the 
situation where p < p < pq, {1~D)L^ ^ 1, and L is much 
larger than the correlation length for given p and D. As 
shown before, the system with (1 — D)L^ <C 1 is well de- 
scribed by the MF model. Since L is assumed very large 
and this value of p corresponds to the active phase of the 
MF theory, the MF model, and accordingly the TC itself 
with above mentioned parameters, stays at the steady 
state of the MF theory with nonzero density for long time 
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FIG. 4: (color online) Schematic representation of three sce- 
narios regarding the limiting value of pc as D — > 1. po indi- 
cated by the open circle is the mean field transition point and 
scales in this figure are arbitrary. The broken lines depict the 
anticipated phase boundaries of each scenarios with known 
transition points depicted by the solid curve. The limiting 
value p of each scenario is the point where the broken line 
meets the vertical line with po on it. The first scenario (SI) 
is not acceptable because the steady state density should in- 
crease as D increases (see the text). The second scenario (S2) 
which was proposed in Ref. is possible only when there 
is a line of critical points (see the text), which does not seem 
plausible. Hence the last scenario (S3) is concluded to be the 
right one. 



of the order of 0(exp(L)) (note that time in this context 
is the rescaled time t). On the other hand, the system 
size is assumed much larger than the correlation length 
and the system is in the absorbing phase by assumption, 
so the density should decay exponentially after time of 
the order of 0{L'^), which is contradictory to the previ- 
ous consideration. Of course, if the correlation length in 
this regime diverges faster than (1 — D)^^^^ as D ^ 1, 
the assumption (1 — D)L^ <C 1 with L larger than the 
correlation length is not valid. However, this scenario of 
the diverging correlation length in this regime also does 
not seem plausible, because it suggests that the corre- 
lation length should diverge for p < p < po and should 
become finite as soon as p > pq. Hence, we conclude that 
p is equal to po- The above discussion is summarized in 
Fig.H 

Up to now, we have not resorted to any assumptions 
about the transition nature. If we assume that the tran- 
sition is always continuous, what can we say about the 
limiting behavior of the TC at p — p = po- In this 
case, the correlation length diverges as (l — D)^'^^ , where 
iy± ~ 1.09 is the critical exponent of the DP class. Ac- 
cordingly, it is not possible to think about the system 
size which is much larger than the correlation length but 
(1 — D)L^ ^ 1 when D is very close to 1. That is, the 
fast diffusion and thermodynamic limit do not commute 
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FIG. 5: (color online) Plot of pc as a function of D for the TC. 
The MF transition point po is indicated by an arrow. Inset: 
Plots of ln(po — Pc) vs ln(l — D) and its fitting function (see 
text). Symbols are from the simulations and the straight line 
is from the fitting. 



at p — pq. However, even if we assume that there is a dis- 
continuous transition for finite D, the noncommutability 
of two limits at p — po is still applicable, because the 
steady state density for p — po is zero for any value of 
D < 1, though the MF theory has a finite density at 
p = Pq. In a sense, this noncommutability is originated 
from the fact that the MF theory exhibits the discontinu- 
ous transition which is characterized by the discontinuity 
of the density at stationarity. At any rate, the thermo- 
dynamic limit commutes with the fast diffusion limit in 
most cases except at p = po- 

Although the above consideration reveals that the 
phase boundary should approach to the MF transition 
point, it does not give any information about the transi- 
tion nature. To get a nontrivial conclusion, the numeri- 
cally obtained phase boundary will be examined. 

Figure O depicts the phase boundary in D — p plane. 
As argued before, the transition point increases with D. 
In the range < D < 0.8 where no controversy has 
ever been raised, the critical points do not change much, 
compared to the change from D = 0.8 to D ^ 0.98. It 
is likely that the phase boundary approaches to the MF 
transition point with infinite slope. If we fit the phase 
boundary using the fitting function 
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Hpo -Pc) = - ln(l -D) + b, 



(3) 



with two fitting parameters 4> and 6, we get Pc ~ Po — 
0.057(1 - i:))°-i27 ^ g-) jj^g^ ^-jjj.gg points [without 

Po and with pc — 0.08630(1) for D — 0.9; see Inset of 
Fig. [5] . Although the accuracy of the fitting should not 
be exaggerated, the clean power-law behavior shown in 
Inset of Fig. [5] strongly suggests that the phase bound- 
ary approaches to the MF transition point with infinite 
slope. Note that this infinite slope is the characteristics of 
the crossover behavior [g^] • Also note that a "crossover" 
within a single universality class, which is actually not a 
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crossover at all, does not show such a singular behav- 
ior [l^. Although we are considering the discontinu- 
ous transition, it is a natural generalization of the claim 
in Ref. fl9'| that no singularity can appear in the phase 
boundary along which only discontinuous transitions oc- 
cur. Thus, the singular behavior of the phase boundary 
near D = 1 supports that no discontinuous transition 
exists for D < 1. 



V. SUMMARY AND DISCUSSION 

In summary, we investigated the triplet-creation model 
(TC) in one dimension. By extensive numerical simula- 
tions, we refuted the previous estimation of the tricritical 
point [ll|, . We only observed the directed percolation 
scaling up to D = 0.98 which is larger than the tricrit- 
ical point suggested in Ref. [Ill, [l^l • To go further be- 
yond D > 0.98, we analyzed the phase boundary near 
D = 1. At first, we argued that the phase boundary 
should approach to the mean field transition point as 
D 1. Using this information, the phase boundary was 
analyzed using the power law fitting function Eq. ^ to 
find ~ Sjwhich, according to the general theory of the 
crossover [20| , strongly suggests the absence of discontin- 
uous transitions in one dimension. 

As a final remark, we would like to comment on the 



conserved ensemble of the TC. In Ref. [ll|, Fiore and 
de Oliveira studied the conserved ensemble of the TC. 
They convincingly argued that the conserved version is 
equivalent to the TC studied in this paper, based on the 
proof in Refs. [U, [22. This argument also embraces the 
limiting case {D 1); see Sec. Ill in Ref. [11]. Ac- 
tually, the equivalence becomes trivially true for the MF 
model under the thermodynamic limit because the equiv- 
alence criterion which is Eq. (30) of Ref. [23| becomes 
identical to the steady state condition of Eq. ([T]). Thus, 
the direct coinparison of the results in this paper with 
those in Ref. [ll| is fully legitimate although we studied 
the different ensemble. In this respect, this work clearly 
suggests a caveat. Unlike the belief regarding the merit 
of the conserved ensemble (for example, see Sec. V in 
Ref. [IH), one should study the conserved version with 
caution even if the discontinuous transition is the main 
interest. 
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